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^v^j We derive the rate of change of the mean motion up to the second post-Newtonian 

order for inspiralling compact binaries with spin, mass quadrupole and magnetic dipolc 
^-j i moments on eccentric orbits. We give this result in terms of orbital elements. We also 

present the related orbital phase for circular orbits. 
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Observations by Earth-based gravitational wave observatories are under way aiming 
to detect gravitational radiation. Upper limits from interferometer data were already 
(f-) . set on inspiral event rates for both binary neutron stars 1 and binaries of 3 — 20 solar 

mass black holes. 2 The parameters of spinning compact binaries can be estimated 
and alternative theories of gravity can also be tested from these measurements. 3 

An important characteristic of these binaries is the rate of decrease of the orbital 
period T due to the energy and angular momentum carried away by gravitational 
waves. Here we give the radiative change of the mean motion n = 2n/T (for eccentric 
orbits). We also present the related change occured in the orbital phase (for circular 
orbits). In both expressions we include all known linear perturbations for an isolated 
compact binary. These are the post-Newtonian (PN), spin-orbit (SO), spin-spin 
(SS), self spin (Self, quadratic in the single spins), quadrupole- monopole (QM) and 
■ magnetic dipole-magnetic dipolc (DD) contributions. 

The expression of the radial period, defined as half of the time elapsed between 
the turning points, emerges from generic considerations on the perturbed Keplerian 
motion. 4 ' 5 Collecting all linear contributions the mean motion has the following form 



Gm 



where rj = fi/m is the ratio of the reduced mass \i to the total mass m of the binary 
system, and £ — —Ej\i where E is the conserved energy. Remarkably there are 
no explicit spin, quadrupolar and magnetic dipolar contributions in the functional 
form of the mean motion. These however contribute implicitly to n through £. 
Since the mean motion is a function of E alone, its evolution can be computed as 
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(dn/dt) = —l/n(dn/d£)(dE/dt). All linear contributions to the secular energy loss 
(dE/dt) due to finite size effects are explicitly given in the literature, 6-9 in terms of 
dynamical constants. The PN contribution is also well-known. 10 Employing these 
we find the change of the mean motion: 

dn\ , T 

dn\ Gm r . T , T , 

^L=- 56cMl- e 2 ) ^ +N >^> < 3 > 
dn\ G 1 / 2 

It) 80 = "APF^ [N * Sl + N ^ ' (4) 

{t) Se = ^a>ll-e*Y^ ,(5) 

/dn\ S^z 
\ dt /SS ™32c 2 m^a 2 (f - e 2 ) 2 

x [Ni sin/ti sinK 2 cos2(i/> — VO + N% cosm cosk 2 + Ng cos 7] , (6) 

dn\ m 2 2 

7t7/ = T~2?~i 2\2^2 Pi [^io(2-3sin 2 K 4 ) + A^iisin 2 K 4 cos2(-0 o -'0 4 )] ,(7) 

at / QM 4a V 1 e J i=1 

$) DD °- ^t^ W ° A + y ' A> ' (8) 

where we have introduced the notations Sl = (S1COSK1 + S2COSK2), = 
[(m 2 /mi) S\ cos /ti + (1111/1112) S2 COSK2] and 

Here Kj and are the polar and azimuthal angles of the i th spin vector and the 
numeric coefficients 
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The orbital elements a, e were derived 11 from the turning points of the ra- 
dial part of the perturbed motion cf. r max = a (lie). (In these variables n — 

(Gm/a 3 ) 1/2 [1 + (77 - 9)Gm/2ac 2 ] ). 

For Kcplerian motions the orbital frequency lo — n. Due to the perturbations, 
precessions occur in the plane of motion (like periastron advance), and the plane 
of motion can also evolve. Therefore the relation of uo and n becomes more com- 
plicated. 12 In the presence of the PN, SO, SS, Self, QM and DD perturbations, 
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for circular orbits (e = 0) the change in the orbital frequency due to gravitational 
radiation is: 13 ' 14 



duj\ 96r/m 5 / 3 w n / 3 
~dt 
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, (10) 



where ct = cr ss + cr Se;/ + ctqm + &dd and /?, fr Sl s 2 , cse//, vqm, &dd arc the spin- 
orbit, spin-spin, self-interaction spin, quadrupole-monopole and magnetic dipolc- 
dipole parameters. 14 For completeness we have added the 2PN and tail contribu- 
tions 15 and we note that higher order contributions are also known. 16 ' 17 

In terms of the dimcnsionless time variable r = r\(t c — t)/5m, defined in terms 
of the time (t c — t) left until the final coalescence, the accumulated orbital phase 
is 4> = 4>c — (5t7j/?7) J w(t)cIt, where <j> c is an integration constant. To second post- 
Newtonian order: 
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The modification induced by the finite size effects SS, Self, QM and DD are all 
encoded in a, while the SO contribution is in (3. 
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